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On Long Waves, 

By J. H. M. Wedderbuen. 



The present paper deals with the theory of long waves in a canal of vari- 
able section. As it is easily shown* that any case of variable section can be 
replaced by one in which the depth alone varies, it is assumed throughout that 
the breadth is constant. 

The main object of the paper is to obtain special cases in which the hydro- 
dynamical equations admit of simple solution for the case of progressive waves ; 
for most of the known soluble cases, apart from the case of constant depth, refer 
either to motions of an oscillatory character or to a reflected train of waves. 
It is also shown that the problem of waves in a basin which is symmetrical 
about a vertical axis reduces in a simple manner to the one-dimensional case. 

The ordinary approximate dynamical equations are assumed without any 
attempt to discuss their validity, although, so far as I know, this has never been 
done in a completely satisfactory manner. From the mathematical standpoint, 
then, our problem is merely the discussion of the solution of a certain partial 
differential equation. The treatment of the initial conditions which is given 
involves the discussion of a type of expansion in series which is in itself of 
considerable interest and seems capable of further extension. 

1. The Equations of Motion. Let the #-axis be taken along the canal 
and the y-axis vertically upwards, the origin being at the undisturbed level of 
the liquid. As usual the displacements, parallel to the axes, of a point on the 
surface are denoted by £ and tj, while h = h(x) represents the depth reckoned 
positive for points below the undisturbed level of the surface. 

For long waves the dynamical equation then isf 

and the equation of continuity is 

8(fc£) dw , 

ox ox 

where 

w=ft£. (3) 

* See for instance Lamb, " Hydrodynamics," 3d ed., p. 257 ; and Chrystal, " On the Hydrodynamical 
Theory of Seiches," Transactions of the Royal Society of Edinburgh, Vol. XLI, p. 614. 
t Cf . Lamb, loc. cit., p. 240. 
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212 Weddeebubn: On Long Waves. 

If vi is eliminated between (1) and (2), we get 

d 2 w _ 3 2 w 
W~ 9 dtf' 



^^; ( 4 ) 



and if £ is eliminated, 

dt 
which, if we set 7* = \dx/h, becomes 

3 2 j? g d' 2 v/ 



?='»{*D' (5) 



a* 8- ^ 8 ' (6) 

which is of the same form as (4). 

Although equations (4) and (5) have both been frequently used, the paral- 
lelism between them seems generally to have escaped notice. It is, however, 
of importance; for if w is a solution of (4) for a certain law of depth h=%(oc), 
then w x (X,) = — dw/dx is, when expressed in terms of X, a solution of 

d' 2 w 1 _ . d 2 w 1 

where ^(X) = l/%(x) ; so that w i (x)is a solution when the law of depth is 
h=4>(x). A repetition of this process evidently leads to the original solution, 
so that the relation between w and w 1 is a reciprocal one. They will therefore 
be called reciprocal solutions. 

2. Wave Propagation in Two Dimensions. Let the x- and 0-axes be taken 
on the undisturbed surface, and the #-axis vertically upwards as before. If a, 
b and c are the initial coordinates of the point x, y, z, and £, vj, £ the corre- 
sponding displacements, the Lagrangian equation of continuity, 

d{x,y,z)/d{a, b, e)=l, 
becomes for long waves 

db ~ \Ba "*" dcJ' 
Since £ and £ are supposed independent of the depth, we have, on integrating, 

yi= ^ {b+h) (|| + <|y_ Ma+ £ c+ £) +Ma> c) , 

which, if a = hl;, y = h%, gives to a first approximation 

\3a 3c/ 3a 3c 

To the same order of approximation the Lagrangian equations for a point 
on the surface are 
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&£_ _ _ dr, d% _ d n 

dt 2 - g da' dt 2 ~ 9 dc 
or 

d 2 a 7 %n d 2 y dy , Q > 

w=- gh r a > w=- gh dc- (8) 

Eliminating yj by means of (7), 

a^^Ha^ + ai-J' # = Ha^+a#> (9) 

or eliminating a and y, 

&='(£*! +£*£)• (l0 » 

If the motion is symmetrical about the ?/-axis, we may set 

» a d c 

r r 



where R is a function of r = Va 2 + c 2 alone. On making this substitution we find 

d 2 R , /3 2 .R . 1 dR 



(11) 



3 2 -R_ /d 2 R ldR ± \ 
af 2_fi, Ur 2 + rar r^J' 

a 2 >7 _ / a , a>7 ^ a>7\ 
W~ g \d~r <k + 7d~r)' 

which may also be put in the form 

d*rR . ,.&rR d 2 n d f A , i dn\ , 19 v 

which are of the same analytical form as the equations for one-dimensional 
waves. This latter form of the equations is derived directly if we consider 
the basin replaced by an indefinite number of canals each lying between two 
vertical planes which pass through the central axis, the angle between the 
planes being infinitesimal. 

Besides the reciprocal transformation suggested by (6), the equations (12) 

are susceptible of another transformation which leaves the form unchanged. 

Let s = l/r; then 

d 2 R 

dt* 

This suggests the more general homographic transformation discussed in the 

next paragraph. 

3. Generalizing the transformation indicated above by combining with it 

a change of scale and origin, let us set 

ax+b , 1 

x' = '—-,, w'= r-j w '> 

cx + d' cx + d 



. /d 2 R , ldR 1 _\ . jd 2 R , 1 dR 1 _\ 
= gh {w + rW-T 2R ) =ghs \W + lW-¥ R )' 
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then, if A=ad — be, 








d Ad 


3 2 


A 2 3 2 


2cA a 


dx~~ {cx+d) 2 dx" 


3a; 2 


(cx+d)*dx' 2 


(cx+d) s dx' 


and therefore 









d 2 w 3 2 dw' d 2 w' 

W 2= W (cx+d)w' = 2c^ + (cx+d) -^r 

2cA dw' 2cA dw' A 2 d 2 w' 

7 i 



{cx+d) 2 dx' (cx+d) 2 dx' ^ (cx+d) s dx' 2 



so that the equation becomes 



where 



aw ,,d 2 w' , 1Q . 

~dF =h -dx^' (13) 

h'= {cX '~ a)i h, (14) 

h being expressed in terms of x'. The transformed equation (13) being of the 
same form as the original equation, we can derive by this method the solution 
for a considerable range of forms of canal from one particular case. 

4. Certain Soluble Cases. Generalizing the solution for the case of con- 
stant depth, we shall now inquire under what conditions (4) has a solution of 

the form 

w=f(x)F[$(x)±ot], 

F being an arbitrary function, but / and definite functions of their arguments. 
Substituting this value of w in (4) we get 

ghf"F+gh(2f'$'+fd")F'+f(ghO' 2 —o 2 )F"=0, (15) 

whence, if F is arbitrary, 

/"=0, 2f'$'+fd"=0, gh$' 2 —a 2 =0, (15') 

which on being integrated give 

f=ax+b, 6=l/a(ax+b), h=o 2 (ax+b) i /g, (16) 

where a is supposed not zero and one unnecessary constant of integration is 
suppressed. 

When a — 1 and &=0, the corresponding values of £ and ri are 



yi=f(- ±oi)—-F'(- ±ot). 
\x / x \x / 



(17) 
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The reciprocal solution is more interesting, as it gives a law of depth more 
nearly realized in practice. Here 

7, = $dx/h=g$dx/o 2 x 4 = — g/3o 2 a?, 
4,{-K)=g/oW=C6 i oW/ g y, 
Q{x)=l/x= — {Zo 2 X/g)K 

Hence, with a slight change of notation, we have the solution 

g=x-*F(at±at)—x- 1 F'(rf±at), 

when h = 9o 2 x*/g. 

The form of (18) suggests assuming a solution of the for 

w=2f n F™(d±ot). (19) 

As before, this leads to the following equations : 

fo=0, 
f" + 2f o 0'+f o d"=0, 



(20) 



/; +2/;_ 1 0'+/„_ 1 0" =f.-,(£-v)> 



Here 6 may obviously be chosen arbitrarily ; but, if the series (19) is to ter- 
minate, we must have 6' 2 = a 2 /gh. 

We shall now investigate some interesting particular cases in which the 
series does terminate. 

The first equation of (20) gives f =ax+b. There are two cases according 
as a is or is not zero ; and it is interesting to find that these two cases are 
reciprocal. 

Let us first suppose that a=f=0. Then, by a simultaneous change of origin 
and scale, we may set f =x. If now we assume 6=x n , it is easily found that 

f r =k r x m+1 , 

Jc r being calculated by means of the recursion formula 

r(rn + l)k r = — [(2r—l)n+l]k r _ 1 , k Q = l. 

In order that f s+1 may be zero, we must have 

= 2f s 6'+f s 6" =k s [2(sn+l)n+n(n— l)W+»*-\ 

which is satisfied if n= — l/(2s + l). This gives 

h=— (2s + l)V», (21) 



216 



Wedderbubn: On Long Waves. 



where w=4(s + l)/(2s + l). The following table gives the values of n and m 
for 5 = 1, 2, . . . ., 5. 



s 


m 


— n 





4 


1 


1 


8 


i 


2 


12 


i 


3 


1 6 

~7~ 


i 

T 


4 


20 
~9~ 


1 
•J 


5 


24 
TT 


1 
TT 



If a=0, we may assume / =1, when we readily find that f r =k t x Tn ; and the 
series terminates after 5 terms if n = 1/(25 + 1). This gives 



h = — (2s + l) V, 
9 



m — 



4s 



25 + 1 



(22) 



5 


m 


« 








1 


1 


4 


1 


2 


8 
7 


1 
T 


3 


1 2 


1 


"7" 


T 


4 


1 6 


1 
9 


5 


TT 


1 
TT 



In the first case the series of values of m approach the limit 2 from above, and 
in the second case they approach the same value from below. For m—2, it is 
easily seen that the two reciprocal solutions coincide. 

By the transformations given in §3, we can derive the solution from these 
results when 

h=(ax + by(cx + d) m /(ex+f) m , (23) 

the coefficients being arbitrary and m having the same values as above ; and also, 
as a particular case, the solutions in the first case can be derived from those of 
the second case and vice versa. 

5. The Normal Form of the Differential Equation. Though by no means 
necessary, it is convenient to bring equation (4) to its normal form. 
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Let d=$dx/Vgh, and put 

r=6+t, 8=6— t, k(r+s)=—6"/±6' 2 ; 
then 

d z w , / dw , dw ' 



=<£+£)• <24 > 



3r3s \ 3r 3s 

It is interesting to note that this equation is the same as that used by Eiemann 
in the theory of waves in air of finite amplitude. 

To obtain a solution of this equation we assume 

CO 

w= 2 *»„[*•<»> (r) +#<">(«)], (25) 

n=0 

where F(r) and G(s) are arbitrary functions of r and s respectively, and m is 
a function of r+s. Substituting this value of w in (24), we get 

»o' — 2km' —0, 



m'^—2km' n =~m' n _ 1 +km n _ 1 , 

whence, if we set m=2m n z n , m satisfies the differential equation 

m"+(z— 2k)m'— kzm=0, (26) 

where the independent variable is r-\-s and z is an arbitrary parameter. 

In discussing the initial conditions in § 9 below, we shall find that it is 
convenient to replace r+s by a variable q=(r+s)/2. We therefore set* 
m(r + s)=v(q), whereupon the differential equation becomes 

£ + „(._,, *_2«=0. (27) 

where 

x(q)=2k(r+s). 

6. A Study of v as a Function of z. In discussing the solution of equa- 
tion (27) we shall assume that q is a real variable lying in a finite range 
So<S<2o+C> an d that x is a real function of q alone which remains finite and 
in absolute value less than K in the range considered, while z is an arbitrary 
parameter which may be imaginary. Under these conditions it is always pos- 
sible to expand v as a power-series in z, and we shall now show that for certain 
determinations of the initial conditions this series represents an integral func- 
tion of z. 

Let v=2v n z n ; then 

v':-2xv' =0, v';-2xv' n ^-2(v' n _ 1 -xv n _ 1 ) ; (28) 

♦When it is desired to emphasize the fact that v depends on a, it will be denoted by v(q, z). 
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therefore 

, \2Kdq C {lain 7 

v =e ) , v o=J e dq, 
<=-2e^ , /fl-^V<_ 1 - W ._0dfl=-2fl^ , Je-^^(«-K._ 1 )d a . 
By a particular choice of the constants of integration we have 

v' =be^ q « iKd9 =ba, v Q =$ q v' dq+a, 
where a=e'*> * Q , and 

v' 1 = -2af^ Tq (v /Va)dq+c, i\ = $ qo v[dq, (29) 

< = — 2a J* y^^K-i/Va)^, v .= J»< rf 2> 

where the constants are so chosen that v n =0, (%>1), when q = q . 
Let now a denote a mean among the values of a ; then 

Vov 
and similarly for n>l, 



v[ = -2^+2a^=r+c = 2Ji(-v +aVa)+c, (30) 






(30') 



where of course the value of a depends on n. 

Let iV and N be the smallest and largest values of \v' \, and set 
L — 2'\N /N ; and let N n and M, denote the maximum values of \v' n \ and \v n \ 
respectively. It must be noted that L is independent of b, so that 2 \ a/a < L. 
Using this notation, we readily derive from (29) and (30) that 

2ST 1 <L(jJf +|fla*| + -^-), ^i<£(^o+|aa*l + -^-)k-2o|, 

#„ < LM n _ x , M n < LM n _ x | q-q 1 , 



whence 
where 

Therefore 
which gives 



M=(M +|aa*|++^-). 
Jf„<L»Jf|g-g |»/», 
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which in turn leads to 

ff„<LMf|<z-gol"-7(»-l): 

M n <L n M\q—q \ n /n(n—l), 
and finally 

N n <L n M\q— g |" _ 7(n— 1)!, M„<L n |2— g | ■/»!. 
We therefore have 

\v\ <M2L n \q— q \ n \z\ n /n\=Me L ^- qM . (31) 

Hence Me L ^~ q "^ is a dominant function for v which therefore converges for all 
values of z, and consequently represents an integral function of that variable. 
The choice of initial conditions made above gives v = a, v' — b + cz when 
q=q . The reason for this choice appears when we consider the equation 
when x is constant. The solution in this case is 



v— e<"-*>*[^e9 v * i! +'< 2 +#e-« Ve2 +'< !i ], 

where for the sake of simplicity q =0. The simplest way in which this value 
of v can be made uniform is by putting A=B, which leads to the same kind of 
initial conditions as the above. There are, of course, other ways in which v may 
be made uniform. 

For reasons which will appear later we make c= — a, so that, when q = q , 

v' =b, v =a, v[=—a, v 1 =0, v' n =0 = v n , (»>1), 

v' — b— az, v=a. (32) 

7. The Asymptotic Expression for v. The asymptotic expression of the 
solution of a differental equation, regarded as a function of a parameter, has 
been the subject of several important memoirs;* but the form in which the 
results are obtained does not appear to enable one to identify the various possible 
asymptotic forms with solutions which are defined with reference to initial con- 
ditions of the type here employed. The method used below avoids this difficulty, 
but it is not applicable to the most general form of equation. 

Let p = qz; then equation (27) becomes 

d 2 v -A x\dv 2xv 
dp 2 \ z/ dp 

For sufficiently large values of z, v can be expanded in the form 



+<!-*)£--=«■ (33) 



f(p,*)=i%, (34) 

n=.0 Z 

where 

<+2vJ=0, v: + 2v' n =2x(v' n _ 1 +v n _ 1 ), (35) 



* Cf ., for instance, Birkhoff, Transactions of the American Mathematical Society, Vol. IX (1908) , p. 219, 
where full references are given. 
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the dashes now indicating differentiation with regard to p. Solving these 
equations, we get 

v' =—2ae- 2 P, v =ae- 2 ?+(3, 

a and (5 being arbitrary constants ; and also 

v n -2Se~ 2 ^e 2p x{v' n _ 1 +v n _ l )dpdp. 

We shall first consider the equation which results from putting for x the 
maximum value, K, of \x\. The solution of this equation is 

V(p, z )=e^- 1 )^-^[Ae^-^ T¥ ^^+Be-^-^ T + w ^}. (36) 



If A+B and (A — B)/\l-\-K 2 /z 2 are uniform integral functions of 1/s, then V 
is also uniform and integral, so that we may write V='ZVJz n , and we can de- 
termine the limits of integration so that V is derived from 

v' — — 2ae- 2(2, - p °>, v =ae- 2< - p - p "^ + b, 

v' 1 =ce- 2 ^-^ + 2e~^ m e 2 *>x(v'o+v )dp, v^i^vidp, (37) 

v' n =2e- 2 ^ Pni e 2p x(v' n _ 1 +v n _ 1 )dp, v n =S PM v'Jp 

by putting K for x. 

If the constants of integration are real, then V is real when p and s are 
real ; we shall assume that p nl and p n2 do not increase without limit as n in- 
creases, so that we can always choose a real value of p so large that V n and 
V'„ have the same sign as V' n _ 1 -\-V n _ 1 . For instance, we can put A and B equal 
to the same real constant a, which gives 

V(p ,z)=2a, V'(p ,z)=(~-l)2a; (38) 

while, if p n2 =p Q (»>1), P nl =Po (»>1), we get from (37) 

V(p ,s)=a + b, 

V'(p ,z) = -2a+ C -, 

which is identical with (38) if o—a and c is chosen equal to 2aK. If we 
therefore give c this value, the function V, defined by (37) when x=K, is 
identical with the function given by (36). Now 

V^b + ae-^-^, V'o+V^b—ae- 2 ^-^ ; 
therefore, if a and b are positive and a<b,p>p and p>p nr , (r= 1,2), for every 
n, then cV and V' -\-V o are positive, and hence V' n and V n are positive for every 
value of n, except w = 0, when V' is negative. 
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Returning to the general case, we shall suppose the limits of integration 
and the arbitrary constants identical with those just chosen for V. The range 
of values of p with which we are concerned is the range of qz, q being real and 
lying between fixed limits, while z is constant but so large as to make \p\ 
greater than p nl and p n2 for all values of n. Let p = qz=re ie ; then 

K| <2Ke- ircqs ^ rni e 2reose \v' n _ 1 +v n ^ 1 \dr. 
Assume now that 

| <_!+«._! | <U!Ci(r cos 6) +V n _ l (r cos 6)], (n>l), 

? ft _! being a positive constant, (Z >1), and — - < 0< - . It follows that 

Therefore l<T*.l<^[r;(f cos c)+F„(r cos «)] 

if r cos is so large as to make V'„, and F n positive when n > ; we therefore 
set / M =^_ 1 /cos 2 6, which gives 

L= k 



cos 2n e 

To determine Z we have 

| v + v' 1 = | b — ae- 2(p - p »> | < Z (6— ae-'C-^ cos *) ; 
if b>a we may therefore put l = (a-\-b)/{a — b), and if a = b a short calcu- 
lation shows that it is sufficient to put Z =l + 2 |tan0|. This leads to 

\v\ <2\v n \/\e\ n <l£r n (r gob 6)/\z cos 2 6\ n <l V(r cos 0, |scos 2 0|), (40) 

since the series for V converges for all values of 1/z different from zero. 
Similarly, for v' we have 

|v'| <Z cos QZV' n {r cos 0)/|ecos 2 0|*<Z o cos V'(r cos 0, |scos 2 0|). 

Now, if z=pe ie , 

■F(r cos \z cos 2 01 ) =ze (K/p cos2 e ~ 1)p(q ~ 9o) cos e [ Ae p(9 ~ Vo) eos ^^^^^^ 

_|_2?g-p(S--4o) cos«Vl+XVP a cos*#i 

p cos 2 0V'(r cos 0, 1 3 cos 2 0| ) = 

e (X/pcos2»-l)p(g-g )cos9r /^" . CQg 2 g i y„2 c S 4 + K Z ) A 6 P (9_9 » ) cos e Vl+tf 2 /P a cos* 

+ (#— p cos 2 0— Vp 2 cos 4 0+^)5e-" (s - 9 »> oostfVl +' K2/ ' ,2cos4( ']. 
These are respectively asymptotic to 

^[gtf(«— 4o)/cos0 an( j J^j^ e K(g—g<,)/coBe 

and therefore | v | and | dv/dq | = | zdv/dp | remain finite as s approaches infinity 

7t It 

in the region 8 defined by — - < < - . 
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When q = q , F(r cos 0, \z cos 2 6\ ) is asymptotic to A + B; and if, in addi- 
tion, A=B, then V (r cos 6, \z cos 2 0| ) is asymptotic to 2AK/p cos 2 6—2A. 

From (40), Z o F(rcos0, s cos 2 0) is a dominant function for v(p, z) as 
regards z, and Z cos 6V (r cos 0, s cos 2 0) is a dominant function for u ' (p, z) , so 



that 



3^ 

dz 

d 2 v 

dpdz 

d 2 v 

dz 2 



<l \dV(r cos 0, p cos 2 0) /dp\, 

<l cos 0|3 2 F(r cos 0, p cos 2 $)/dpdp\, 

<l \d 2 V(r cos 0, p cos 2 0)/8p 2 | ; 



and from these three relations we find, after a short calculation, that zdv/dp, 
z 2 d 2 v/dpdz and z 2 d 2 v/dz 2 remain finite as p approaches infinity. 

If we now assume that A = B = a = b, we have, when p = p or q — q , 



v{p<>,z)=2a, 

But dv/dp = dv/zdq, and therefore 

dv 



dv(p , z) _ c 

dp z 



dq 



= c—2az, (q = q ). 



These initial conditions are of the same form as in (32), b being replaced by 
c and a by 2a. It follows that v is an integral function of z when p is replaced 
by qz. 

Now v, regarded as a function of z, remains finite in the region S defined 

7t 71 

by — - < < „ , and approaches a definite limit as z approaches infinity ; hence 
in the same region zdv/dz and z 2 d 2 v/dz 2 remain less than some finite quantity. 



But 



d 2 v 
dz 7 



d 2 v 



+ 2q 



d 2 v 



d 2 v 
~dz 2 



' 3~2 » 



dp 2 dpdz 

and in this equation z 2 d 2 v/dz 2 , qz 2 d 2 v/dpdz and z 2 d 2 v/dz 2 have been shown to be 
less than some finite quantity in the region 8; hence the same is also true of 
z 2 d 2 v/dp 2 , and therefore of d 2 v/dq 2 . Now 

d 2 v dv fdv \ 

and we may therefore set 



dv _/d 2 v 

dq \dq 2 



)/*=■■-, 



dq)' z 

where U remains finite ass approaches infinity in S. Hence 



Z Jq„ 



Udq. 
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In this expression we know that U is finite in 8 even when q = q ; hence v is 
asymptotic to aeU<, Kdq 

in* S. 

8. On an Expansion of an Arbitrary Function in Terms of v. We shall 
now regard v(q, z) as a function of q and z alone, choosing the arbitrary con- 
stant so that 

., dv , 

dq 
when q = q , so that in 8 

v(q, z)^eU Kdq . 
Now it is readily shown that 

satisfies the same differential equation as v. Differentiating, we have 
u ' = e -^(Q-Qo)[ v '(q } — z )—2sv{q, —z)], 

so that when q = q 

u—1, u' — b + z — 2s = 6 — z. 
Hence us==v, or 

v (Q, -~ z ) = e^o-^v (q, z). (41 ) 

The rest of this paragraph follows closely the treatment of one of Cauchy's 
proofs of Fourier's development, which is given in Picard's Traite d' Analyse, 
t., II, p. 180, to which we refer the reader for some of the details of the argu- 
ment. The method depends on the integral 






i[3f(re*)— 3f(— re*)]e*<ty, (42) 

which is equal to the sum of the residues of the analytic function %(z) inside 
the circle with center at s=0 and radius r. A special form is chosen for $(z), 
and then the limit of (42) is investigated as r approaches infinity over a 
sequence of values so chosen that no circle passes through a singularity of the 
integrand. When this limit exists, it is equal to the sum of the residues of $(z) . 

Cauchy sets 

*(«)/(«, g ) 

7l(z) 

where •$, n and / are integral functions of z, f alone depending on the param- 
eter q. The main difference in the investigation given below is that 4* and n 
are also taken to depend on q ; otherwise the conditions are very much the same 
as in Cauchy's development, being varied only by the special form assumed for 4-. 

* That this is also true on the boundaries of 8 is shown by Birkhoff, Transactions of the American 
Mathematical Society, Vol. IX, p. 219. 
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We shall first of all put 

/ being a real integrable function which satisfies Dirichlet's conditions. 
A sufficient condition for the convergence of the integral (42) is that 

i*^yJ]^- fl) /^)^-i*^=0/'e-^-rt/^)^ (43) 

shall approach a definite limit as z increases, so that its real part remains 

positive and its modulus takes on the increasing sequence of values r 1 , r 2 , 

Departing slightly from Cauchy's conditions, we suppose that ^(z)/n{z) 
approaches a definite limit c, while e^i — z)/n(—z), (b = q — q x ), tends in 
general to zero. The second part of (43 ) can then be written 

and therefore approaches zero under the given conditions.* The first part 
approaches! 

\c lim z f e-°«-> i) f(ix)di*=icf(q). 

where \ is a root of 7i(z), and the summation extends over all the roots. 
Similarly, if b = q 2 — q and 

Jq 

and C is put for c, the same conditions lead to 

itf/(ff)=24^T r%-«/(^)^, (45) 

n (\) Jq 

where of course 4 1 and it are not necessarily the same as before. 
Let us first set 

^(z)=e e< - a+q -^v(q, z)=e^ a - q+q "' ) v(q, —z), 
7 t(s)=e^*-«» ) (e as — 1), 



Then, in 8, 



provided 



v ^(z) n . e»< B+ «-*>t;(2> s) r ,«,„ 
c= hm -4-4 = lim .,„_„>, ng , / =eJg» ,H " > 

lim ^ ( ~ g) «*- lim g^ZT^M^g) _ Q 
!™7t(-s) e - ™ e-<«-*)( 6 — -1) ~°' 

2(S— « )+6<a. 



*Cf. Picard, loc: oit., p. 182-183. fCf. Pieard, Joe. c**., p. 183-184. 
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Hence, if \=2itin/a, (n=0, ±1, ±2, ), denote the roots n{z), we have 

ieU t «*f(q) =\l v(q, \) fW<-»f(fi)dii. (46) 

a _ x J Ql 

Using now the second form of f(z, q), we shall set 

4,(s) =e z( - q -^v{q, z) =e~ z ^-^v(q, —z), 

7t (s)=e 2 («-«°)(l— e~ m ). 



Then 

as before, and 

provided 
Hence 



(7= lim^4=eW dff = c 

n=x 7t(Z) 

lim *$=*>■*•= lim ^'"^ (g ' «> =0, 

2(^— g )+&<a. 

ieW*«/(g)=- 25 «(«,**) f e x » ( « -»f(fi). (47) 

Adding (46) and (47), we arrive at the final form 

eU^f(q) = - 2 v(g, \) f e^-^f((i)d(i, (48) 

provided 

2(3— 2o)+2— 2i<», 

2(2— 2o) + 22—2 <«• 

The form of these inequalities shows that we may put q = without loss 
of generality ; and having done so, we easily find that, if q is any number in the 
interval q 1 <q<q 2 , then the maximum value of q 2 is (2a — qJ/3, while q x <a/2. 
We may, for instance, put q 1 =0, which gives q 2 =2a/3. 

It should be noticed that if q 1 and q 2 lie between and a, then the series 
given in (48) still converges if the factor v(q, \) is omitted from each term, 
as is easily seen by comparing (48) with Fourier's expansion. 

9. Initial Conditions. The only type of initial conditions which we shall 
consider is the case where the profile of the wave and the horizontal velocity 
of particles on the surface are given as functions of x when t=0. Expressed 
in terms of w, this means that dw/dx and dw/dt are given as functions of x 
when £=0. 

Now 

dw _/),/3w dw\ dw __ dw dw 

~dx~ \lh + ~^s) , ~dT~!fr~ ~ds~'' 
we may therefore set 
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dw . dw dty(q) 



dr ds "' dq ' 

(r=s = q), (49) 

dw dw . . . 

where ^ is a given function of # and <£> is determined to an additive constant. 
When r=q=s, the form assumed for w in (25) becomes 

w=iv n {Q){F^{q)+G^{q)} 

= 2 «.(«)#<•> (2), (50) 

where H(q)=F(q) +G(q). Inserting this value of w in the first equation of 
(49), we get 

*'= t («:+«._ 1 )ffW(g) 

« = 

♦=lv.(a)H(->(a); (51) 

K = 

and similarly, if JST(a) =F'(q)—G'(q), 

*= I «._ 1 [FW(a)-(?w(a)]= I «.(<z)£-<*>(<z). (52) 

"We shall now proceed to make use of the developments of the preceding para- 
graph in showing that we can in general determine H and K to satisfy the 
initial conditions. 

Assuming that H and K are expansible in a Fourier series 

H(q)=X a n e^, K(q) = I 0„e^, (53) 

— 00 00 

and that these series can be differentiated term by term, we have, from (51), 
$= I v r (q) f a n Ke K «= la,l ^(g)*^ 



rzz:0 »= — oo *i = — oo r — 



= 2a,«(«A)«'"'. (54) 

— oo 

This expansion is of the same form as (48), where, if 

^ I =e-I«.' ,d «^, (55) 

we have 

-«> (56) 

Giving a n this value in (53), the series derived for If satisfies the required con- 
ditions if these same conditions are satisfied by the Fourier series of the given 
function $. 
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Similarly, for K we obtain 

i Q2 

A=-f«-HW^ (57) 

where ^=e-y*-*. (58) 

From this we easily deduce 

^(2)=;|[ff(<Z)+/ 9 tf(/«)^L (59) 

The form of (56) suggests another form for w which makes it unnecessary to 
introduce functions having derivatives of every order. 
Assume 

w=2v n [i(r+s), \] [a„e x » f +^ M J, (60) 

where the summation extends over any sequence of values \ for which the 
series converges together with its first and second derivatives. It is easily 
verified that this series satisfies the differential equation for w, and the discus- 
sion of the initial conditions then proceeds as before. 

10. Some Special Methods. We have seen in the preceding sections how 
to find a solution for any case in the form of an infinite series, but in some 
particular cases there are more convenient methods. 

Let h = kx; the equation for 57 then is 

3 2 >7 __, 3 / dy;\ 

a particular integral of which is* 

vi=AJ {nx i ) cos not, 
where n is arbitrary and a=^gk. This solution may be written 



7T 

y\=B f 2 cos (nx h sin 6) cos notdO 
•Jo 

B — 
— "o" f 2 [cos n(x* sin 6+ at) + cos n (x i sin 6 — ot)]dd. 



Hence we have as a solution 

7T 

y,= (*[F(at sin 6+ot) +F(rf sin 0—ot)]dS, (61) 

Jo 
F being any even function which can be expanded as series of cosines. 

This gives two wave systems travelling in opposite directions. Each part 
is not, however, by itself a solution. For if we put 

* Cf . Lamb, loo. oit., p. 259. 
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IT 

v,= f*F(at sin 6+at)d0, 

Jo 
we get 

This shows, however, that if G is any odd function, 

W 

17= ( T [G(rf sin 0+<rf) — G(»* sin 6+at)]d6 

Jo 
is a solution, and since any function can be expressed as the sum of an even 
and an odd function, we readily derive a solution in the form 

7T 

J7= C T [F(at sin $+ot)+F(—at sin 6+ot)]d6. 

Jo 
In certain cases it is possible to determine F in (61) to satisfy given initial 
conditions. To show this let ^=x sin 2 6; then, if yj is the value of yj for t = 0, 

/>- /•» F(^) 

»=»j['Wih»)*=j[^* 

A solution of this integral equation was obtained by Abel* in the form 



if / 

»/ 



Vlp — X 





when >7 satisfies the following very general conditions : namely, yj (i) is contin- 
uous, (ii) has a finite derivative in the region under consideration, and (iii) 
vanishes for x=0. 

This is a solution for which dvi/dt (and therefore also d%/dt) is zero when 
t=0. 

Instead of forming the reciprocal solution, we shall investigate independ- 
ently under what conditions there is a solution of the form 

IT 

w=f ( T [F($ sin 0+ot)+F(Q sin 6—ot)]d6, 

Jo 

where / and $ are functions of x. Let 

■w 

Wl =f C*F(<}> sin 6+ot)d6; 

Jo 

substituting in (4), we get 

WW WW 

<r 2 / f T F"dd=gh[f" ( T FdB+ (2f<p' +f<}>" —f<p' 2 /<}>) C*F' sin Odd+fQ' 2 C J F"dO] 

Jo Jo Jo Jo 

+ 9 -^fF'{°t), 

__> 

* See, for example, Bdcher, " Introduction of the Study of Integral Equations," 1909, p. 6. 
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whence 



/"=0, 

P 

gh$' 2 —a 2 =0, 



2f'^'+fr-~=0, (62) 

<P 



and F' is an odd function. 

Integrating the first two equations, we have 

f=av+b, j=WS («*0), 

= k, (a=0). 

Taking first the case a =0, let /=1; then $ = e kx , the constant of integration 
being absorbed in a, and 

a 2 

Zi — p —ilex 

h ~ gk 2& ' 

£=^—e 2kx C T [F(e k * sin $+ot) +F(e kx sin 0— at)]dd, (63) 

w Jo 

TT 

V = —Jce kx C T [F'(e k * sin 6+at) +F'(e kx sin 0— at)] sin 0d!0. 

This solution is the reciprocal of (61) except that k must be replaced by ^k. 
If a:£0, we mayHby a change of scale and origin set f=x, which gives 

q> = e- k/x , h=a 2 x i e 2k/x /gk 2 . 

We have already found that there is a solution w = f (x) F (0 ± at), 
F being arbitrary, only when h= (ax-\-b) i . If, however, F is not arbitrary, a 
solution of this type may be found for any form of canal. 

Since the coefficients in (15) are functions of % alone, while F is a function 
of 6+at, this equation leads to 

ghfO' 2 —a 2 f _ 1 



ghf" a 

and 

gh(2fd'+fd") 



(64) 
= /?, (65) 



ghpd' 2 — <T 2 / 
a and /3 being constants. F must therefore satisfy the equation 

F"+pF' + aF=0. 
We shall first consider the case /?=0. In this case we have 

F = C cos (i(d±at), (i 2 =a; 
also 

2/'0'+/0"=O, 
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whence 

B'f = h, 
and, from (64), 

/At 2 P 
g (t 2 k 2 — f f 

so that when h is given / is determined. In practice this equation for / is 
intractable, but interesting cases arise on giving / particular forms. For 
instance, if f=x n , we get 

(l 2 2 X in 

h= 



and for n=\ 

0= log x. 

The latter case is interesting as giving a particular solution of the self- 
reciprocal case previously referred to. 

Another soluble case is obtained by setting /= (% 2 +bx+c) n . 

If instead of putting /?=0 we assume gh$' 2 — <j 2 =0, or l/a = 0, the equation 
to be satisfied by F becomes 

" P+ 2f6'+f6" F=0, 
and 

/" 



2fd , +f6 r ' = constant =y 

or 

/"— 2y/'0'— y/0"=O, 

the normal form of which is 

V 2 (7 2 

v"+^-v=0, f=e«v. 
gh 



Hence 



w=ve ± v". 



which is equivalent to the ordinary harmonic solution of (4). The equation 
can also be normalized by changing the independent variable x to x 1 — ^e iye dx, 
which gives the solution 

w=v(^ $ dx)e- y( - e± '' t \ 



Princeton University, June, 1912. 



